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SUMMARY 

This  document  discusses  the  use  of  the  Discrete  Fourier  Transform  to  estimate  the  phase  of 
particular  frequency  components  of  a  sampled  signal.  Attention  is  paid  to  the  limitations  of 
the  Discrete  Fourier  Transform  in  general.  Particular  emphasis  is  placed  on  estimating 
phase  for  use  in  Interferometric  Direction  Finding. 
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1  INTRODUCTION 

The  principle  of  Interferometric  Direction  Finding  is  based  on  the  difference  in  phase  of  a 
signal  frequency,  measured  at  two  physically  separated  antennas.  This  paper  explores  the 
possibility  of  using  a  Discrete  Fourier  Transform  (DFT)  to  obtain  phase  from  a  signal,  and 
hence  the  prospect  of  using  the  DFT  In  Direction  Finding  (DF)  applications. 


2  INTERFEROMETRIC  DIRECTION  FINDING 

The  bearing  of  a  radio  emitter  may  be  determined  by  phase  interferometry  (1).  This  involves 
the  use  of  two  or  more  spatially  separated  antennas  in  an  array.  The  phase  difference 
between  the  signals  received  at  each  antenna  can  be  used  to  determine  the  angle  of  arrival. 

Consider,  as  an  example,  the  two  clement  interferometer  shown  in  Figure  1. 


Figure  1  Two  element  interferometer 

Due  to  the  difference  in  path  length,  the  phase  difference  of  the  signal  on  arrival  at  each 
antenna  is 


=  22E_l  =3  L  = 

X  2* 

The  bearing  of  the  emitter,  9,  may  then  be  computed  as 
0  =  arccos  (L/D) 

NOTE  This  is  a  simple  arrangement  and  will  result  in  an  ambiguity  tn  the  direction  of 
arrival,  about  the  antenna  baseline.  In  practice,  at  least  three  elements  placed  in  a  triangle, 
or  four  in  a  crossed-baseline  arrangement  are  required  to  remove  the  ambiguity. 
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The  phase  difference  between  the  signals  received  at  each  antenna  must  also  be  constant  in 
time.  Thus,  the  frequency  of  the  component  received  at  each  antenna  must  be  exactly  the 
same.  In  practice,  this  means  that  when  a  receiver  employs  a  frequency  conversion  stage  prior 
to  sampling,  the  same  local  oscillator  must  be  used  for  all  antenna  charnels. 


3  THE  DISCRETE  FOURIER  TRANSFORM 


The  Fourier  Transform,  X(f)  of  a  continuous  signal,  x(t)  is  defined  as 


X(f)=  x(t>e-l2Kft.dt 


Ihc  analogous  DFT  of  the  signal,  sampled  at  N  points,  T  seconds  apart  is  given  by 


N-l 


X(k)  =  £  x(n)e-'2'dtn/N 

n*0 


N-l 

=  £  x(n)  (cos(2rtkn/N)  -  jsin(2rckn/N)l 

n»0 


where  k  and  n  are  the  discrete  equivalents  to  f  and  t,  related  by  t  =  nT  and  f  -  k/NT. 

NOTE  Both  n  and  N  arc  required  to  be  integers,  while  k,t,T  and  f  are  real  numbers. 

X(k)  is  a  sequence  of  complex  numbers  giving  the  discrete  spectrum  of  x(n).  The  quantisation 
steps  in  the  spectrum  will  be  referred  to  as  frequency  "bins".  From  these  complex  numbers,  the 
magnitude  and  phase  of  each  point  in  the  spectrum  may  be  obtained  using  simple 
trigonometry. 

The  Fast  Fourier  Transform  (FFT)  is  simply  a  method  of  computing  the  DFT  with  much  fewer 
calculations  than  the  direct  method,  and  is  thus  much  faster  to  compute.  Most  algorithms  for 
the  FFT  require  that  the  number  ot  points,  N,  be  restricted  to  a  power  of  two. 

3.1  Relative  phase  in  the  Fourier  Transform 

As  phase  is  a  relative  quantity,  it  is  important  to  know  what  the  phase  of  the 
transform  X(f)  is  measured  with  respect  to. 

The  phase  of  a  particular  frequency  component  of  a  waveform,  as  determined  using  the 
Fourier  Transform  is  measured  relative  to  the  first  sampling  point,  with  the  equation 
describing  the  waveform  being  a  COSINE  relationship,  as  shown  in  Figure  2. 


2 


UNCLASSIFIED 


UNCLASSIFIED 


ERL-0530-RN 


Figure  2  An  example  of  the  relativity  of  phase  of  a  sinusoid 

As  an  example,  the  following  functions  have  "phases"  resulting  from  the  Fourier 
Transform  as  follows: 

SIN  (Ft +  0)  •-=>  DFT  phase  =  -90°  ,atf  =  F 

SIN  (Ot  +  90)  =>  DFT  phase  *  0°  (  DC  signal ) 

SIN  (Ft  +  90)  =»  DFT  phase  =  0*  ,  at  f  =  F 

3.2  The  need  for  windowing 

The  DFT  operates  on  a  finite  length  data  stream,  or  sample  period.  If  the  frequency  of  a 
particular  component  of  a  waveform  is  such  that  an  integral  number  of  period?  of  that 
component  is  contained  in  the  sample  period,  then  that  component  has  no  associated 
frequency  error.  This  means  that  when  the  DFT  of  the  waveform  is  taken,  the  frequency 
of  that  component  falls  exactly  into  one  frequency  bin.  If  however,  the  period  does  not 
have  this  property,  then  the  component  will  not  fall  exactly  into  one  bin,  ar.d  that 
component  is  said  to  have  a  frequency  error. 

The  frequency  error  of  a  component  is  the  amount  by  which  the  frequency  of  that 
component  deviates  from  the  frequency  of  the  closest  bin.  All  real-world  signals 
exhibit  frequency  error. 

The  effects  of  the  frequency  error  on  the  magnitude  of  the  spectrum  have  been  well 
documented  [2,31  and  arc  generally  known  as  Spectral  Leakage  and  the  Picket  Fence 
Effect.  Spectral  Leakage  occurs  when  the  frequency  of  a  component  does  not  fall  exactly 
into  one  bin,  and  refers  to  the  spreading  of  the  energy  from  one  frequency  bin  to  adjacent 
bins. 

The  Picket  Fence  Effect  draws  an  analogy  between  the  output  of  the  DFT  of  a  waveform, 
and  a  bank  of  analogue  bandpass  filters,  each  filter  being  centred  at  one  of  the  frequency 
bins  in  the  discrete  spectrum.  Ideally,  each  "filter"  should  have  a  square  frequency 
response,  but  in  practice,  the  responses  are  rounded,  and  the  result  is  that  the  spectrum 
appears  to  bo  viewed  through  a  picket  fence. 
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The  Spectral  Leakage  and  Picket  Fence  effects  are  worst  when  the  frequency  lies 
midway  between  two  frequency  bins.  ie.  a  frequency  error  of  0  5, 

In  order  to  reduce  these  effects,  "windows"  are  applied  to  the  sampled  data  prior  to 
performing  the  DFT.  Windowing  is  the  multiplication  of  the  waveform  by  a  discrete 
weighting  function  -  the  window. 

3.3  The  effect  of  windowing 

Consider  a  signal  x(t),  windowed  by  the  window  function  w(t),  to  produce  the  signal 
y(t).  The  signal  x(t)  has  a  Fourier  Transform  X(f>,  and  w(t)  has  the  Fourier  Transform 
W(f). 

y(t)  =  x(t)  x  w(t) 

Then,  Y(f)  =  X(f)  *  W(f)  C  denoting  convolution) 


|X(X)  WU-iOM^^I.dX. 


In  the  discrete  case,  this  integral  becomes  a  summation. 

N-l 

Y(k)  =  £  |X(X) 

As  ar.  example,  consider  a  cosinusoidal  signal,  x(t>,  sampled  at  N  points  and  having 
some  discrete  time  window  w(n)  applied  to  it. 

x(t)  a  Acos(2jtf0  +  4>) 

it  can  be  shown  that  the  continuous  Fourier  Transform  of  x(t)  is  given  by 

^  X(f)  =  A.  5(f-f0)  ci*  +  A.  §(f+f0)  e-i* 

2  2 

where  5(f)  is  the  delta  function. 
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It  follows  that  the  DPT  of  x(n)  is  then 

X(k)  =  A.g(k-ko)  d*  +  A-^k+ko)  e‘l* 

2  2 

The  DFT,  W(k)  of  w(n)  can  be  expressed  in  the  form 

W(k)  =|W(k)|cie<w 
N-1 

V (k)  =  £  X  (k-X)  W(X) 

X-0 

N-1 

=  £  A.(8(k-X-ko)ei*  +  «k-X+ko)c-l*>|Wtt)|d*« 
o  2 

N-1  .  .  N-1 

=  X  eK*9(M1£-5<k-X-k0)|  W(X)!+  X  d('**4KU>  A.  8(k-X+ko) |  W(X)| 

o  2  o  2 

But,  5(k-X-k0)  =  0  at  all  X,  except  at  k-Xko=0,  when  5  =  1,  and  8<k-X+k0)  =  0  at  all  X, 
except  at  k-X+ko=0,  when  5=1. 

V(k)  =  A.|  w(k-k0)|  oi(*  ♦“‘-W  +  A I  W(k+k0)| 

2  2 

3,3.1  The  boundary  effect 

Equation  2  consists  of  two  terms  which  arise  as  a  result  of  the  window  function 
acting  on  the  positive  and  negative  frequency  components  of  X(k) 
respectively. 

it  is  customary  to  simplify  the  equation  by  ignoring  the  effect  of  the  second 
term.  However,  blind  use  of  this  simplification  may  lead  to  incorrect  results, 
as  the  contribution  of  the  second  term  to  positive  frequency  points  is  non-zero 
in  the  general  case.  The  real  part  of  the  positive  frequency  spectrum  appears 
to  be  reflected  about  the  zero  frequency  point,  while  the  imaginary  part  of 
the  spectrum  is  reflected  with  sign  inversion.  For  a  finite  length  sample 
sequence,  this  reflection  also  occurs  about  the  point  N/2.  Hence,  the  points 
close  to  N/2  are  also  similarly  affected. 

The  extent  of  the  contribution  of  negative  frequencies  will  be  largest  for 
signals  with  a  frequency  close  to  0  or  fs/2  (k=N/2),  windowed  by  functions 
with  large  sidclobes  in  the  frequency  domain. 

Figure  3  shows  the  linearity  of  the  magnitude  from  a  128  point  FFT  of  a 
discrete-frequency  signal  with  a  frequency  error  of  0.4.  Each  point  corresponds 
to  the  magnitude  of  the  output  of  an  FFT  performed  on  a  sinusoid  at  that 
frequency.  A  rectangular  time  window,  which  has  a  sinc-like  Fourier 
Transform  with  large  sidclobes  was  used  throughout.  It  can  be  seen  that  the 
lowest  and  uppermost  12  discrete  frequencies  suffer  errors  of  more  than  1  %  in 
magnitude.  A  similar  sized  phase  error  is  also  incurred. 
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Figure  4  shows  the  linearity  improvement  when  using  a  Hanning  time 
window.  This  window  has  a  much  narrower  Fourier  Transform  and  as  such 
only  the  first  and  last  two  points  suffer  errors  above  1  %. 

The  effect  may  thus  be  reduced  by  applying  suitable  windows  and/or 
considering  points  far  away  from  the  spectral  extremes.  Using  longer  length 
transforms  allows  a  greater  number  of  these  points  to  be  considered.  The  effect 
may  be  eliminated  by  applying  the  Hilbert  Transform  to  the  sampled  data  io 
obtain  the  "analytic"  signal.  Such  a  signal  contains  only  positive  frequency 
components. 

For  simplicity,  the  remainder  of  the  paper  will  assume  negative  frequency 
components  have  negligible  effect,  unless  otherwise  stated.  Under  such 
conditions,  equation  2  simplifies  to 

Y(k)  A.|W(k  -ku)|  . .,(3) 

From  equation  3,  it  can  be  seen  that  for  a  single  frequency  waveform,  the 
magnitude  of  the  spectrum  will  consist  of  an  amplitude  scaled  Fourier 
Transform  of  the  window  function,  centred  at  the  frequency  of  the  signal 
component  k^. 

If  the  signal  has  no  frequency  error,  then  the  spectrum  of  the  window  will  be 
centred  exactly  on  a  frequency  bin.  Thus  the  best  window  to  use  in  this 
circumstance  is  the  rectangular  time  window.  As  this  has  a  Fourier  Transform 
of  a  sinc-like  function,  there  will  be  nulls  at  all  bins  other  than  at  it’s  centre. 
The  transform  then  will  consist  only  of  one  line,  at  the  frequency  of  the  signal. 

By  considering  equation  3  again,  it  can  be  seen  that  the  phase  of  the  spectrum 
will  be  the  phase  of  the  Fourier  Transform  of  the  window,  centred  at  the 
signal  frequency,  plus  the  phase  of  the  signal  component. 

Again,  for  a  signal  with  no  frequency  error,  using  a  window  with  nulls  at  most 
frequency  bins  other  than  it's  centre  will  result  in  the  phase  from  the  DFT 
being  a  good  estimate  of  the  true  phase  of  the  spectrum.  This  is  because  if  the 
leaked  amplitude  response  is  very  small  at  a  particular  bin,  then  the 
contribution  to  the  overall  phase  at  that  bin  wilt  be  small. 

NOTE  If  the  magnitude  is  zero  (or  very  small )  at  a  frequency  bin,  then  phase 
does  not  exisi  for  that  point  and  can  be  arbitrarily  set.  in  these 
investigations,  the  phase  in  this  case  was  set  to  zero. 
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3.3.2  Multiple  frequency  components 

Most:  typical  signals  used  in  DF  are  real  signals,  containing  no  imaginary 
component,  They  may  therefore  be  treated  as  a  set  of  frequency  components, 
each  having  particular  amplitude  and  phase. 

For  the  case  where  the  signal  x(t)  comprises  more  than  one  frequency 
component,  the  resultant  Fourier  Transform  will  be  a  superposition  of  the 
Fourier  Transforms  of  the  constituent  components. 

ic.  Y(k)  *  £  Ajw(k-k()1)!expj((j>i  +  6(k-ko())  . (4) 

l  2 

Although  it  would  be  tempting  to  think  so,  this  does  not  mean  that  the 
resultant  phase  and  magnitudes  are  each  themselves  superpositions  of  the 
constituent  phases  and  magnitudes. 

The  real  case,  with  no  frequency  error  will  have  the  Fourier  Transform, 
evaluated  at  some  particular  frequency  kn,  given  by 

Y(kn)  =  A  ,/  +  O(0))  +  A  V  |w(k„-ki)|  Z  (pi  +  9(kn-ki>)  . (5) 

2  2  i X 

If  the  window  used  is  rectangular,  then  equation  (5)  simplifies  to  the  expected 

result  Y(kn)  <=  A  z<t>n. 

2 

This  indicates  that  with  a  rectangular  window,  and  no  frequency  error,  each 
component  has  a  distinct  phase,  and  no  component  affects  any  other 
component. 

3.3.3  Experimental  justification 

The  result  for  magnitude  is  expected,  and  is  well  known.  The  phase,  on  the 
other  hand,  seems  intuitively  satisfying,  but  it  must  be  shown  that  this  result 
agrees  with  results  from  actual  FFT’s.  Fourier  Transforms  have  been 
performed  on  128  point  sampled  sinewaves  at  frequency  f0-32,  and  using  three 
windows  -  Rectangular,  Triangular,  and  Hanning,  with  frequency  responses  as 
follows: 

Rectangular 


W(k)  =  c->™^— 1^-L 

jsin^yNj 
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Triangular 


Hanning 


n  ,  sinM 

W(k)  =  - 1-2-L 

N  sln(itk/N> 


w(n)  = 

W(k)  =  c'I^'t1’  (0.5D(k)  +  0.25  [D(k  -  1/N)  +  D(k  +  1  /N>]  > 


where  D(k)  =  e'*lKk/N  ■■ _ — 
sin{nk/N) 


The  magnitude  and  phase  resulting  from  performing  an  FFT  on  the  sampled 
waveform  using  each  of  the  above  windows,  has  been  plotted  and  appears  in 
Figures  5  to  10.  The  theoretical  phase,  as  calculated  using  equation  (3)  is 
marked  with  an  "X". 


NOTE  There  is  only  one  non-zero  point  in  the  spectra  produced  using  the 
Rectangular  window,  and  three  points  using  the  Hanning  window  -  ihe  centre 
‘oh:,  and  the  t wo  adjacent  bins.  I'nis  is  because  nulls  jrom  the  frequency 
response  of  the  windows  fall  in  all  other  frequency  bins. 


0.00  6.40  12.8  19.2  25.6  32.0  38.4  44,8  51.2  57.6  64.0 

FREQUENCY 

Figure  5  Magnitude  of  FFT  of  SIN,  f  »  32,  using  a  rectangular  window 

The  Rectangular  and  Hanning  cases  have  shown  exact  correlation  with  the 
theoretical  values. 
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FREQUENCY 

Figure  9  Magnitude  of  FFT  of  SIN,  f  3  32,  ^  using  a  triangular  window 
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The  Triangular  ease  has  shown  very  good  correlation  for  frequencies  near  f0, 
but  accuracy  falls  off  as  frequency  moves  further  from  f0.  This  is  to  be 
expected,  as  the  magnitude,  at  distant  frequencies  is  very  small,  and  thus  a 
ratio  of  two  very  small  numbers  is  taken  to  obtain  phase.  Any  error  in  either 
the  real  or  imaginary  component  will  cause  large  phase  error. 

The  Rectangular  and  Hanning  windows  do  not  suffer  from  the  same  errors  as 
the  Triangular  window  because  the  centres  of  the  frequency  bins  fall  at  the 
nulls  of  the  magnitude  spectra  of  these  windows. 

3.3.4  The  effect  of  windowing  on  the  phase  from  the  Discrete 
Fourier  Transform 

Considering  equation  (5),  it  can  be  seen  that,  for  the  purely  real  case  with  no 
frequency  error,  the  phase  at  kn  is  dependant  on  the  following  factors: 

1  the  phase  of  the  component  at  frequency  kn. 

2  the  phase  of  all  other  frequency  components  kj. 

3  the  phase  at  fn  of  the  window  function,  centred  at  each  of  all  other 
frequency  components  kj. 
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NOTE  The  extent  to  which  factors  2  and  3  above  affect  the  phase  at  kn,  is 
proportional  to  the  magnitude  of  the  window  function  centred  at  k{,  at  kn-  To 
reduce  the  effects  of  factors  2  and  3,  the  window  function  should  be  chosen  such 
that  W(k)  is  minimal  at  frequency  kn-  When  there  is  no  frequency  error  the 
optimal  window  to  use  is  the  Rectangular  window ,  because  the  Fourier 
Transform  of  this  window  has  nulls  at  all  frequency  bins  ether  than  it's  centre. 

3.3.5  The  effect  of  frequency  error  on  the  phase  from  the  discrete 
Fourier  Transform 

Most  typical  signals  have  components  with  frequency  error.  In  this  case, 
equation  (2)  still  holds,  with  kn  and  k|  being  the  exact  frequencies  of  the 
components  -  net  the  frequency  bins  in  which  they  fall. 

This  has  a  detrimental  effect  on  the  use  of  most  window  functions.  For 
example,  consider  the  rectangular  window.  In  the  case  of  no  frequency  error  a 
rectangular  window  has  spectral  nulls  at  all  bins  but  the  centre.  However, 
with  frequency  error,  the  spectrum  of  the  window  will  no  longer  be  centred 
exactly  at  one  bin.  The  result  is  that  the  nulls  no  longer  full  at  bin  centres  • 
there  will  now  be  finite  magnitude  levels  at  the  bins.  Hence  there  will  also 
be  a  phase  contribution  at  most  bins. 

As  stated  in  paragraph  3.3.4,  the  effects  of  the  contribution  by  window 
functions  to  magnitude  and  phase  can  be  minimised  by  reducing  the  magnitude 
of  the  Fourier  Transform  of  the  window  function  at  the  appropriate  frequency 
bins.  This  implies  the  use  of  windows  with  low  sidelobe  levels. 

An  example  of  the  effect  of  frequency  error  on  the  phase  of  the  spectrum 
resulting  from  a  DFT  is  shown  in  Figures  11  to  14.  Using  a  sampling  frequency 
of  128  Hz,  a  128  point  transform  has  been  performed  on  a  sinusoid  of  frequency 
f0  =  32  Hz.  The  sinusoid  has  been  windowed  by  using  a  triangular  window. 

Figure  11  shows  the  phase  of  the  spectrum  with  no  frequency  error.  At  k=32, 
the  phase  can  be  seen  to  be  the  correct  -90° .  Figure  12  shows  the  phase  of  the 
spectrum  with  the  sinusoid  having  a  frequency  error  of  0.1  (ie.  ko  =  32.1).  The 
phase  at  k  -  32  is  seen  to  be  approximately  -72°. 

Figures  13  and  14  show  that  the  phase  at  k  =  32  is  reduced  to  -48®  and  0®  for 
frequency  errors  of  0.2  and  0.5  respectively.  Thus,  it  can  be  seen  that  the 
phase  at  32  becomes  progressively  less  correct  as  the  frequency  error  increases. 
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Figure  12  Phase  of  1  FT  of  SIN,  f  =  32,  with  frequency  error  =  0.1,  using  a  triangular  window 
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FREQUENCY 

1  igure  13  Phase  of  1  IT  of  SIN,  f  =  32,  with  frequency  error  =  0.2,  using  a  triangular  winuow 
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Figure  14  Phase  of  FFT  of  SIN,  f  -  32,  with  frequency  error  =  0.5,  using  a  triangular  window 
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4  DIRECTION  FINDING  AND  THE  DISCRETE  FOURIER  TRANSFORM 

When  using  the  DFT  for  DF  applications  the  difference  in  phase  between  a  frequency 
component  of  a  signal  received  at  one  antenna  and  that  received  at  another  antenna  must  be 
obtained  with  high  accuracy. 

The  DF  calculation  using  physically  separated  antennas  may  be  simulated  by  performing  a 
DFT  on  a  waveform,  altering  the  phase  of  the  frequency  components  of  that  waveform  (to 
simulate  wave  propagation)  and  then  performing  another  DFT  on  the  new  wave. 

4.1  Single  frequency  component  signals 

Examination  of  equation  (3)  reveals  that,  tor  a  single  component  signal,  a  DF  type 
phase  difference  calculation  should  yield  exact  results,  as  the  window  dependant  term 
will  cancel  in  subtraction. 

To  verify  this  observation.,  a  simulation  was  performed  which  involved  constructing  a 
4096  point  discrete  sinewavc  at  frequency  Icq,  zero  phase,  and  then  performing  an  FFT  on 
it.  The  phase  of  the  sinewavc  was  then  altered  by  90°,  and  the  FFT  was  repeated.  The 
phase  deference  from  each  FFT  at  the  frequency  k0  was  obtained,  and  the  percentage 
error  calculated.  The  simulation  was  performed  for  all  values  of  lq,  between  0  and  2048, 
and  the  percentage  error  plotted  against  frequency  as  shown  m  Figure  15. 


FREQUENCY 

Figure  15  Linearity  over  frequency  of  phase  difference  error 
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A  rectangular  time  window  was  used  with  each  frequency  ko  having  a  frequency  error  of 
0.4,  as  this  gave  the  most  conservative  results.  Figure  15,  shows  that  the  error  is 
constant  and  near  zero  at  all  points  except  those  close  to  the  ends  of  the  spectrum.  The 
deviation  at  these  points  is  due  to  the  Boundary  Effect,  and  may  be  minimised  as 
explained  in  paragraph  3.3.1. 

The  nonlinearity  towards  the  edges  of  the  spectrum  means  that  the  bandwidth  over 
which  accurate  relative  phase  estimation  may  be  obtained  is  slightly  less  than  the 
expected  bandwidth  of  0  to  fs/2.  If  multiple  bands  are  to  be  used  to  analyse  a  wide 
bandwidth,  then  some  overlapping  of  these  bands  may  be  necessary  to  obtain  the 
desired  linearity. 

A  second  simulation  was  performed  to  verify  the  linearity  over  actual  phase  difference 
of  the  calculated  phase  difference.  Some  128  point,  32  Hz  sinewaves  were  generated  at 
integral  phase  differences  from  0  to  360°.  FFTs  were  performed  on  the  sinewaves,  and 
the  calculated  phase  differences  were  plotted  against  actual  phase  differences,  refer 
Figure  16.  Again,  a  rectangular  time  window  was  used,  and  the  sinewaves  had  a 
frequency  error  of  0.4,  to  allow  the  most  conservative  results  to  be  obtained.  The  plot 
shows  almost  perfee*  linearity. 


FREQUENCY 


Figure  16  Linearity  over  actual  phase  difference  of  calculated  phase  difference 
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4.2  Multi-component  signals 

Referring  back  to  equation  (4),  one  could  draw  a  similar  conclusion  to  that  of 
paragraph  4.1,  that  the  window  dependant  term  will  cancel  on  subtraction  in  a  DF  typo 
phase  calculation  of  a  multi-component  signal.  However,  this  assumption  is  incorrect. 
The  window  dependant  term  of  equation  (4)  is  also  dependant  on  <j>j,  the  phase  of  the 
ith  frequency  component  of  the  signal.  The  parameter  <j>i  will  not  be  the  same  when 
measured  at  spatially  separated  antennas,  due  to  propagation  of  the  signal.  Thus,  the 
window  dependant  term  does  not  cancel,  and  so,  the  resulting  phase  difference 
calculation  is  not  exact. 

The  use  of  windows  with  low  sidelobcs  reduces  the  error  due  to  the  phase  contribution 
by  other  "distant"  frequency  components,  kj.  However,  most  common  windows  such  as 
the  Hanning  window,  have  a  main  lobe  in  the  frequency  domain  that  is  wider  than  one 
frequency  bin,  and  thus  any  adjacent  or  nearby  frequency  components  will  still  contribute 
to  the  phase  at  the  frequency  of  interest,  kn.  This  is  tolerable,  provided  the 
amplitudes  of  these  components  arc  small.  However,  if  the  amplitude  of  adjacent  or 
nearby  components  is  of  the  same  order  as  the  amplitude  of  kn,  the  phase  corruption  is 
too  great  to  allow  any  reasonable  accuracy  to  be  obtained  with  a  DF  type  phase 
calculation 

Thus,  simple  direction  finding  using  the  DFT  is  possible  only  for  narrowband  signals. 
The  tolerable  bandwidth  is  dependant  on  qualities  of  the  window  function  such  as  main 
lobe  width  and  sidelobe  amplitude  levels. 


5  RESULTS 

The  following  experimental  observations  have  been  made  on  the  phase  obtained  using  the 

DFT  to  process  purely  real  signals.  These  observations  are  qualified  by  equations  (2)  and  (4). 

It  is  expected  that  these  observations  should  apply  also  to  the  case  of  complex  signals  (ie. 

those  possessing  both  real  and  imaginary  parts). 

1  A  DF  type  relative  phase  measurement  of  a  single  frequency  component  signal,  with 
or  without  frequency  error  will  achieve  correct  results,  provided  that  allowance  be 
made  for  the  Boundary  Effect  -  the  contribution  of  negative  frequency  components  to 
the  positive  part  of  the  spectrum. 

2  A  DF  type  relative  phase  measurement  of  a  multiple  frequency  component  signal, 
with  frequency  errors,  will  not  achieve  exact  phase  estimation.  Phase  corruption  by 
other  frequency  components,  as  well  as  that  due  to  the  negative  part  of  the  spectrum 
occurs.  However,  the  use  of  windows  with  low  magnitude  sidelobcs,  and  narrow  main 
lobes  will  ensure  a  reasonably  accurate  estimation. 
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6  CONCLUSION 

The  Discrete  Fourier  Transform  is  capable  of  providing  an  accurate  estimate  of  relative 
phase,  as  required  for  its  use  in  direction  finding.  Windowing  sampled  data  prior  to 
performing  the  DFT  reduces  both  the  Boundary  Effect  and  the  effects  of  frequency  error  on 
relative  phase  measurements,  and  thus  greatly  enhances  accuracy.  Results  of  the 
investigation  have  indicated  that  the  use  of  the  Discrete  Fourier  Transform  for  direction 
finding  is  feasible. 
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